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Abstract 

We give a general survey of the solution of the Einstein constraints 
by the conformal method on n dimensional compact manifolds. We 
prove some new results about solutions with low regularity (solutions 
in H2 when n=3), and solutions with unsealed sources. 

Dedicated to Vincent Moncrief, a great scientist and friend. 



1 Introduction. 

The geometric (physical) initial data on an n-dimensional manifold M for the 
Einstein equations are g a properly riemannian metric and K a symmetric 
2-tensor. They cannot be arbitrary, they must satisfy the constraints, which 
are the Gauss-Codazzi equations linking the metric g induced on M by the 
space time metric g, and the extrinsic curvature i^' of M as submanifold 
imbedded in the space time {V,g), with the value on M of the Ricci tensor 
of g. 

These constraints read as equations on M. They are the so called hamil- 
tonian constraint (we choose units such that StcG = 1): 

R{g) - K.K +{trKf = 2p, (1.1) 

and the momentum constraint: 

V.K - VtrK = J. (1.2) 
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p is a scalar and J a vector on M determined by the stress energy tensor of 
the sources and the lapse N . In a Cauchy adapted frame, where the equation 
of M in y is = t =constant, one has: 

with p > if the sources satisfy a positive energy condition. 

In this article we use the standard method^ to solve the constraints, that 
is the conformal method, initiated by Lichnerowicz and developped by CB 
and York. We treat the cosmological case^, that is we suppose the manifold 
M to be compact (without boundary). We consider arbitrary dimensions 
and we lower the regularity previously obtained for the solutions: we obtain 
large classes of solutions in H2 when n — 3. We prove also miscellaneous new 
results. There is certainly still room for improvements. 

2 The conformal method. 

A riemannian metric 7 is arbitrarily chosen on the manifold M, together with 
a scalar function r. The physical metric g and second fundamental form K 
are then defined by^: 

4 . . — 2n(n+2) . . 1 • ■ 

gij = ip^^ij, K'' = ip^^A'' + -g'^T, (2.1) 

with A"^^ a traceless symmetric tensor. The sources are supposed to be known, 
they split into York scaled and unsealed ones, they are given as a pair of non 

^The constraints were first formulated as an elliptic system in CB 1957, explained also in 
CB 1962, through harmonic spacetime coordinates. The decomposition of the constraints 
through the conformal method into a linear system and an elliptic semi linear equation 
was obtained for zero trace K and zero momentum sources by Lichnerowicz 1944. An 
elliptic formulation for the linear system was given in CB 1970. The decomposition of the 
system of constraints was generalized to constant trace K and scaled momentum by York 
1972, who obtained an elliptic system by the now universally used conformal splitting. 
Improved variants of this splitting, useful for numerical computations, have been given by 
York 1999 and Pfcistcr and York 2002. 

^The asymptotically euclidean case (isolated systems) is considered in CB- Isenberg- 
York 2000. Results with similar low regularity can be proved by analogous methods. They 
arc being obtained independently by D. Maxwell who also treats the case of asymptotically 
euclidean manifolds with boundary.. 

•'In the case n=2 the powers of if are to be replaced by exponentials (see Moncrief 
1986). 
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negative scalars Pi, P2, and a pair of vectors Ji, J2- The physical sources are: 

2(n + l) 2(n + l) 

q = ip + J = (p "-2J1 + J2. (2.2) 

One denotes by D, and by A^, the covariant derivative, and the Laplace 
operator, in the metric 7. 

The Hamiltonian constraint becomes the Lichnerowicz equation: 

n = A^ip-f{.,ip)^0, (2.3) 

/•/ \ _3n-2 j!^ n+2 

j[.,ip) = Tip - aLp "-2 -gi</? "-2 + (6 - 0^2) 9? "-2, 
with (a dot is the scalar product in 7) 

'n — 2 ^, . n — 2,, , n — 2, 

r=- -i? 7 a = - -A.A>Q, 6 = -— > 0, 

4(71-1) 4(n-l) 4n 

n-2 

The momentum constraint becomes (indices raised with 7): 

M' = DjA'^ - {^^^^^d'r + J2 + Ji} = 0, (2.4) 

In the original conformal formulation one sets: 

A^^ = B'^ + {C,,^fXf (2.5) 

with B'^^ an arbitrary given traceless tensor, and C^^conf^ the conformal Lie 
derivative of a vector X to be determined, that is: 

{C^^,^fX)ij = DiXj + DjXi - --i,PhX\ (2.6) 

irt 

We will call the system 2.3 and 2.4, like Isenberg and Moncrief, the LCBY 
equations. The unknowns are ^p and X. 

Remark 2.1 In the thin sandwich formulation (York 1999), which makes 
the conformal invariance of the constraints (when r = constant) more trans- 
parent and is useful in numerical computations, one sets: 

A'i = {2N)-\-u'^ + {C,,confPf} (2.7) 
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with N an arbitrary, positive, scalar, and u^^ an arbitrary given traceless 
tensor. The momentum contraint in this formulation reads, for general n: 

M' = Dj{{2N)-\C,,,,^,fPy^} F^si-, ^) = 0, (2.8) 

Fi^si; V) = Dj{{2N)-^u'^} + !^(^2n/n-2^i^ ^ ^2(n+2)/(n-2) ^ 

The unknowns are now (3 and ip. The mathematical results are essentially 
the same as in the original formulation which we use here for simplicity of 
notations. 

Functional spaces. 

To have good functional spaces it is convenient to introduce a given 
smooth properly riemannian metric e on M. 

The Sobolcv spaces on (M, e) are defined as closures of spaces of 
smooth tensor fields in the norm 

•^^ 0<fc<. 

where d, \.\ and /Xg are the covariant derivative, the pointwise norm and the 
volume element in the metric e. A metric 7 is said to belong to if it is 
properly riemannian and 7 e W^. We will always suppose that a > -, then 
7 e is continuous on M as well as its contrevariant associate 7*, and 
MP is an open set in W^. The volume element fj,^ of 7 is equivalent to //g. 

3 Solution of the momentum constraint. 

The momentum constraint 2.4 reads: 

{A^,confXy = D^{C^,,^fXf = F\ (3.1) 

pi = !illl(^2n/(n-2)5V + ^2{n+2)/(n-2) ^ _ .g ^) 

n 

where S is a traceless symmetric 2-tensor. When 7, properly riemannian, r, 
Ji, J2, and B are given and ip is known 3.1 is a linear system for the vector 
X. 



4 



Lemma 3.1 1/7 e Mf^^, s > 0, p > | the vector F belongs to VFf , as soon 
ascpe Dt, JeWP andBe W^+^. 

Proof. The Sobolev multiphcation theorem. ■ 

Theorem 3.2 The momentum constraint 3.1 with 7 e Mf^^, p > |, s > 
/las a solution X e WI+sj ^/ £ ^'^^ -^^ orthogonal to the space of 
conformal Killing (CK) vector fields of (M,j). 

The solution is determined up to addition of a CK vector. It is unique 
if we impose that it be orthogonal to CK vectors. There exists then a 
constant > depending only on 7 such that 

II X ||w|V2< II F ■ (3.3) 

Proof. We show that the operator A^,conf satisfies the theorems 7.2 and 
7.4 of the appendix A. 

1. The operator is eUiptic: its principal symbol at x, with ^ G T^M is the 
linear mapping from covariant vectors b into covariant vectors a given by 

C%bj + eCjbi - hj^bk = % . (3.4) 

This linear mapping is an isomorphism if ^ 7^ : its characteristic deter- 
minant is computed to be (7* is the contrevariant tensor associated with 
7): 

Dio = ictrii - -) = 7*(c, > if e 7^ 0. (3.5) 

n 

2. The coefficients of the operator A^^conf are of the type : 

^2=7*, 01=7*^7, 00=7*9797 + 7*9^7 (3.6) 
they satisfy the hypothesis of the theorem 7.2 of the appendix A if 

02 e M|, ai e WP, ao G with P> ^ ■ (3.7) 

Indeed if 7 G M|,p > f , then also -j* e d C°'°^, while 97^7 e 
Wf X C LP, hence Oi e Wf, and Oq G i/. 

2. Smooth vectorfields X and Y with smooth 7 satisfy, on a compact 
manifold, the following identity: 

/ y,(A^,eon/X)X = / Y.dAd'X^ +D^X'--^'^WuX^\iJi^= (3.8) 
JM JM \ n J 
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Under the hypothesis made on p the operator Ay^conf is a continuous mapping 
W2 ^ ■ hence F.A^ ^on/-^ is a continuous mapping W2 x W2 ■ The 

conformal Kilhng operator is of the form 5+7* 57, it is a continuous mapping 
from into Wf, and by the Sobolev embedding it holds that Wf x Wf C 
as soon as » > a fortiori if p > ^r. 

The continuity of all the considered embeddings permits the passage to 
the limit which proves the identity 3.8 in our low regularity case. This identity 
implies, making X = F, iS.^^conj^ = 0, the following one: 

(£^,eon/X)^,• = A^,' + ^t^X, - = 0. (3.9) 

The identity 3.8 also shows that /\^^conf — k : VFf is injective, 

hence an isomorphism, if A; is a strictly positive number. Its inverse is then a 
compact operator, hence con/ is a Predholm operator, A^^conf^ — F has 
a solution X e Wf iff F is orthogonal to the kernel of the adjoint operator, 
which is con/ itself. The solution is unique in the Banach space of VKf 
vectors orthogonal to CK vectors. 

The proof of higher regularity for more regular coefficients, and of the 
inequality 3.3 is obtained by derivating the equation. ■ 

Remark 3.3 If the Ricci tensor^ of'-y, with components p^j, is negative def- 
inite the manifold (M, 7) admits no conformal Killing fields. Indeed the 
equality A^^conf^ — implies, using the Ricci identity, 

D'DiXj + (1 - -)DjD'Xi + p.iX^ = 0, (3.10) 

which implies on a compact manifold 

f {-D'XWiXj - (1 - -)DjXW'Xi + pjiX'X^pi^ = 0, (3.11) 

from which follows X = if Ricci(j) is negative definite and n>2. 

Remcirk 3.4 Even when X is not unique (i.e. if{M, 7) admits a CK vector 
field) the tensor A is determined uniquely through the formula 2.5. 

^The hypothesis 7 € , p>n/2 implies that the Riemann and Ricci tensors of 7 are 
in LP. 
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Lemma 3.5 When Dt = and J2 = 0, then F is LP' orthogonal to CK 
vector fields if and only if it is so of Ji. 



Proof. If Dt = and J2 = 0, then F reduces to the sum of Ji and 
the divergence of a traceless symmetric tensor. Such divergences are 
orthogonal to CK vector fields. This fact well known in a smoother case 
results from the identity, still valid under our hypotheses 

/ D.B'^Xjfi^ = - [ B'^{L^,,onsX\,ii^. (3.12) 

JM JM 



If Dt ^ or J2 ^ (unsealed momentum sources) there seems to be at 
present no result for manifolds (M, 7) admitting CK vector fields. 

4 Lichnerowicz equation. 

It reads 

A^(p ^ f{.,(p) = r(p - acp "-2 -qi(p "-2 + (6 - g2)</7"-^ (4.1) 

we consider in this section that the coefficients a,qi,b are given non negative 
functions on M. 

4.1 The Yamabe properties. 

A theorem conjectured by Yamabe and proved by Trudinger, Aubin, Schoen 
is interesting to classify solutions obtained by the conformal method. 

definition 4.1 The functional 

defined for every (p ^ Hi C , (p ^ and 7 e Mf p > |, is called the 
Yamabe functional. 
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The functional Jq{ip) admits an infimum for if G VFf (dense in Hi if 
p > |), and (p ^ 0, because it is bounded below as shown by the inequality: 

/ (/c-i|V(^r + i?(7VK>-i / R{i)^'f^,\>-\\R{i)\y\\'P'\\,^ 

JM JM 

(4.3) 

This infimum /j,: 

H= Inf Jj{<f)= Inf Jy{<f), (4.4) 

depends only on the conformal class of 7, it is called the Yamabe invariant. 
The manifolds (M, 7) are split into 3 Yamabe classes according to the sign 

of /jL. 

definition 4.2 The manifold (M, 7) is said to be in the negative Yamabe 
class if n < 0, in the zero Yamabe class if /i — 0, in the positive Yamabe 
class if fj, > 0. 

Renicirk 4.3 It is known (Kazdan and Warner (1975, 1985) that a sufficient 
condition for (M, 7) to be of negative Yamabe class is J^R{'-f)fi^ < 0, and 
that every compact manifold of dimension n > 3 admits a metric in the 
negative Yamabe class. On the other hand, not all manifolds (Lichnerowicz 
1963) can support metrics in the positive or zero Yamabe class. 

Possible topologies of manifolds of various Yamabe type are reviewed in 
Fisher and Moncrief 1994. 

The following theorem conjectured by Yamabe has been proved in an 
increasing number of cases by Trudinger and Aubin and finally completed by 
Schoen, by showing that the infimum is attained by a function 9?^. 

Theorem 4.4 Let M he a compact smooth manifold. Any smooth rieman- 
nian metric j on M is conformal to a metric with constant scalar curvature, 
the sign of which is a conformal invariant. 

Remcirk 4.5 The Yamabe theorem says that the minimum <^m, hence 7m) 
unique if /i < 0, but not necessarily so if fj, > 0. 
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The Yamabe invariant is defined as soon as 7 e Mf, P > f , since then 
R{'j) G L'P, and (/? G VFf C C°. Checking of the different steps of the proof 
shows that the existence of ip^ G , > can be extended with no 
real difficulty to these low regularity metrics in the cases < and /x = 0. 
Difficulties for this generalization appear in the case fi > and specially 
when n — 3. The proof given by Schoen uses indeed the geodesic balls and 
the Green function which require at least a C^'^ metric, implied only by 
7 G VFf+s! p > f , if s > 2. The weaker theorem proved by Yamabe^ extends 
to 7 G Mf , P > f and will be sufficient for our use. 

Theorem 4.6 Let {M, 7) be a compact riemannian manifold with 7 G Mf , 
p>|. Then: 

1. If is in the negative Yamabe class it is conformal to a metric with 
scalar curvature —1. 

2. If J is in the zero Yamabe class it is conformal to a metric with scalar 
curvature 0. 

3. If 7 is in the positive Yamabe class if it is conformal to a metric with 
continuous and positive scalar curvature r > 1. 

4.2 Non existence and uniqueness. 

definition 4.7 Let f > be a function defined almost everywhere on M. 
We say thaf f ^ if there is an open set of M where Inf / > 0. 

Theorem 4.8 (non existence) The Lichnerowicz equation admits no so- 
lution 9? > 0, (/? G W2, p > J, on a compact manifold (M,j), 7 G Mf, 
a,b,q ^ L^ if either: 

1. r < 0, b — q2 < 0, and — r + a + gi — 6 + ^2 ^ 0. 

2. r > 0, a = gi = 0, 6 - 52 > and r + b-q2^0. 

Proof. The integral of A^^(/9„ with respect to the volume element of 7^ 
is equal to zero on a compact manifold M if 7„ is C-^ and is C^, by the 
Stokes formula. We approximate the given 7 and in C C° by such 7^ 
and </?„. Denote by /i^ and fj,^ the volume elements of respectively 7„ and 7, 

tends in to fj,^. The Sobolev embedding and multiplication theorems 

^See for instance Aubin p. 127. 

^If a is continuous the definition coincides with the usual notation. 
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show that Ajip{fj,^)~^ fj,^ — A^^(p{fj,g)~^ fj,^^ tends to zero in (see analogous 
proofs in the appendix), hence the integral of A^Lp on (M, 7) is zero. 

On the other hand, f{x, (p) is integrable on M since Lp & and > 
and its integral cannot vanish under the condition 1 [respectively 2] which 
implies that /(x, ^) <Q [respectively > 0] and is strictly negative [positive] 
on an open set of M. ■ 

The geometrical origin of the Lichnerowicz equation leads to a general 
uniqueness theorem, independant of the sign of r, hence of the Yamabe class 
of 7. 

Theorem 4.9 (uniqueness^ ) The Lichnerowicz equation 4-1 on (M, 7), 7 G 
^2>P > f ' most one positive solution ip G M^f, if a,qi,b € L^, a>0, 

qi > 0, b — q2 > 0, and a+ qi+ b — q2 ^ on M. 

Proof. Suppose it admits two solutions (pi > and (p2 > 0. The 
following identity holds: 

A,4/(„-.,^(^,^J^) -r(^^/('^-^)7)(^i^2^) ^ -(^i^2^)^"-^'^/^"-'V(^t/(--^)7) 

(4.5) 

Since (p^ is a solution of 4.1 we have an equality of the form: 

r(v^f^"-S) - -v^r^"^'^/^"-'^{A,^,-^ir(7)} 

and an analogous equation for r(<^2'^*^"~^^7)- Inserting in the previous equation 
gives the equation: 



^l^M"^^^ - 1) - M(^i^2"' - 1} = (4.6) 



with 



(,n ,o-1^4(n-l)/(n-2) i 
, _ (-3n+2)/(n-2) -(n+2)/(n-2) (y^l'^2 j ^ ^^ 1 

-n/(n-2) -(n+2)/(n-2) , /, v -l{'^1^2^Y^'"'^~'^^ - ^ 

^1^2 - 1 

The fractions with denominator 'Piip2^ — 1 are continuous and positive func- 
tions on M if it is so of the ip's, because of the mean function theorem which 

^First proved (case n=3, smooth coefficients) by Araki, see CB 1962. 
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says that if a and a are two positive numbers there exist a number b such 
that 

^L_zi ^ ab''-\ 6 e [a, 1] if < a < 1, be [1, a] if a > 1. (4.7) 
a — 1 

Therefore A is an integrable function on M, depending on (p-^, (p^j o^-iQi-iQ^ 

and 6, with A > and A ^ on M if a + gi + 6 - > and ^ 0. The 
equation imphes therefore® ^1^2 ^ — 1 = 0. ■ 

4.3 Existence theorem, scaled sources in L°°. 

The results can be obtained, as in the case n = 3, by using the Leray 
Schaudcr degree, CB 1972 (unsealed sources), O'Murchada and York 1974 
(scaled sources), or a constructive method, Isenberg 1987. Both methods use 
sub and supersolutions. The following theorem is an extension of results ob- 
tained before with higher regularity, and for n = 3 (see reviews in CB-York 
1980, Isenberg 1995). In the theorems some coefficients arc supposed to be in 
and not just in L^, because we look for constant sub and supersolutions. 
A new result is an extension to cases where ^ in the cases of negative 
or zero Yamabe class. 

Theorem 4.10 The Lichnerowicz equation with scaled sources on a compact 
n-manifold (M,e) with given riemannian metric 7 G Mf, p > f , and a, b, 
qi e L°°, with unknown the conformal factor ip, admits a solution ip > 0, 
cpeWi if: 

1. fM, 7) is in the positive Yamabe class and Inf{a + qi) > 0. 

M 

2. (M,j) is in the zero Yamabe class, Infr^ > and Inf{a + qi) > 0. 

M M 

3. (M,^) is in the negative Yamabe class and Infr^ > 0. 

M 

Proof. The general existence theorem 8.2 shows that, under the hypoth- 
esis made on the coefficients, the Lichnerowicz equation admits a solution 
(p > 0, (p e W2, if we can find constants i and m such that on M: 

f{;i)<0, f{-,m)>0 with 0<e<m. (4.8) 

We have: 

^See appendix A, lemma 7.6. 
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f{x,y) = y n-2h{x,y), (4.9) 

h{x, y) = b{x)y^^ + r{x)y "-2' — qi{x)y "-2' — a{x) (4-10) 

There exist £ > and m> i satisfying 4.8 if they satisfy the same inequahties 
with h. We have 

2 ~ n ~ 2"+2 2n-2 

/i' = 7:h{x,y)y'^, h{x,y) = 2nby~^ +2{n- l)ry~^ - (n - l)qi. 

" n — z 

(4.11) 

8 " 

/i^ = -{n(n + l)by^ + (n - lfr}y^. (4.12) 

" n — 2 

The function h'y{x, .) has at most one zero Yq(x) > 0; the function h{x, .) 
starts from — (n — 1)^1(2;) < for y = 0. We examine the various Yamabe 
classes. 

• 1. (M, 7) is of positive Yamabe class, we can suppose r continuous 

and r{x) > 1 at every point x G M. Then for any given a; G M the 
function hy{x, .) is a non decreasing function of 7/ > 0, non negative as 
soon as 

2n— 2 1 

y > m,{x), with m(a;) a number such that m(a;) "-2 > -gi(x). 

(4.13) 

The function /i(x, .) of ?/ is non decreasing when y > m{x), it is non 
negative if it is so for y = m{x). A sufficient condition is that, since 
b > 0,: 

4(n-l) 2(n-l) 

m{x) "-2 > qi(x)'m{x) "-2 + a{x) (4-14) 
which is imphed for instance by the sufficient condition 

m(x) > Max{l, a(x) + qi(x)} (4.15) 

Any number m (independent of x) such that 

m > Sup(a + qi), m > 1 (4.16) 

M 

will satisfy the condition 4.15 almost everywhere on M, hence will be a 
supersolution. 
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The function h{x, .) is non positive for some i{x) if 

in 4(n-l) 2(n-l) 

b{x)£{x)'^ + r{x)£{x)~^ - qi{x)e{x)~^ - a{x) < (4.17) 
A sufficient condition is 

4(n-l) 2(n-l) 

i{x) < 1, {b{x) + r{x))e^^ < {qi{x) + a{x))£^^ (4.18) 
a sufficient condition is tfierefore: 

i{x) < 1, e{x) < {b{x) + r{x))-\qi{x) + a{x)) (4.19) 
Tlierefore if Inf{a + gi) > any number £ sucli tliat 

M 

< £ < {5iip(6 + r)}-^{/n/(gi + a)}, £ < 1. (4.20) 

M M 

is a constant subsolution,and it is possible to choose i and m such that 

0<e<m. (4.21) 
• 2. (M, 7) is of zero Yamabe class, we suppose r{'y) — 0. Then: 

4n 2(n-l) 

h{x,y) = b{x)y"-2 — qi{x)y "-2 — a(x) (4.22) 

There exist £ > and m > £ satisfying 4.8 if they satisfy the same 
inequalities with h. We have 

2 ~ n ~ 2n+2 

/i^ = —-Ji{x, y)y^^ , h{x, y) = 2nby^ - (n - l)q^. (4.23) 

4 n+4 

h!^= n{n-r\)\nr" . (4.24) 
There exist constant supersolutions < £ < m if Injia + 51, 6) > and 

M 

Su-pia + qi,b) < +00. They are then chosen such that: 

M 

Inf{a + qi) Sup{a + qi) 

0<£< Min{l, ^ ^ , }, m> Max{l,— . (4.25) 

Supb Injb 

M M 
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3. (M, 7) is of negative Yamabe class, we take r — —1.. 

An 4(n-l) 2(n-l) 

h{x,y) = b{x)y'^ — y "-2 — qi{x)y "-2 — a(x) (4-26) 
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h'y = ^^h{x,y)y'^-^, h{x,y) = 2nby" ' -2(n- l)y "'-2' - (n- l)gi. 

(4.27) 

/i; = ^{n(n + l)6y'^ - (n - ifjy^K (4.28) 

If b{x) > then /i(a;, .) starts from a non positive value, decreases until 
y^Yo(x), 

increases afterwards up to +00, it has therefore one positive zero Y{x) > 
Yq{x). Hence h{x,.) starts from a non positive value, decreases until 
y — Y{x) then increases up to +00. We see that if Inf 6 > on M one 
can always find £ and m satisfying 4.4 by choosing i and m such that: 

Q<i< InfY{x), m > SupY{x). (4.30) 

M M 



Sub and supersolutions do not need to be constants: we will obtain the 
following theorem by looking for non constant sub or supersolutions. It 
replaces inequahties to be satisfed on M by inequahties on an open set of M. 

Corollary 4.11 The theorem 4-^0 holds if 

1. The hypothesis Inf{a+qi) > is replaced by a+qi ^ for the positive 

M 

or zero Yamabe classes. 

2. The hypothesis Infr"^ > is replaced by: 

M 

a. ^ for the zero Yamabe class. 

b. Inf t'^ > on a sufficiently large subset of M for the negative Yamabe 
class. 
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Proof. 1. We show as in Isenberg 1987 that the condition Inf{a + qi) > 

M 

is superfluous for the construction of a non constant , by constructing a 
constant subsolution for a conformaUy transformed equation. Suppose that 

a + qi ^ 0, i.e. there is an open subset U <Z M where Inf{a + qi) > 

u 

0. The complementary set M — U is contained in a proper open subset V 
which we can take with smooth boundary dV. We can construct a metric 7' 
conformal to 7 whose scalar curvature is stricly negative in V, as follows. Set 
7' — ^^/("~^)7, the formula for the conformal transformation of the scalar 
curvature gives: 

r(7') = e-^{-A^e + r{-f)e). (4.31) 

Take for 9 e the unique strictly positive solution in V of the linear 
equation and boundary condition, with k some positive constant, 

A^^ - A;^ = 0, e \dv= 1- (4.32) 

These two equations imply that: 

r(7') = 6'~^(r(7) - A;). 

We take k > Supr{^) we have then r(7') < 0, with Inf\r{'^')\ > 0, in V. 

M V 

We choose on M a smooth function ^ > equal on V to the previously 
determined 9. The function ip' = 9~^if satisfies the following equation, where 

An An 

a' = a9~"^ , q[ = qi9~^^ , hence inf (a' + q[) > 0: 



n+2 



A^/ip' - r(7')(/?' + a'(f' "-2 + q[(f' "-2 - bf'^-^ = 0. 

The number £' > will be a subsohition of this equation if it is such that (we 
impose i' < 1, which is no restriction, to simplify the writing of the second 
inequality) : 

Inf\r{i)\ Inf{a' + q[) 

0<l'< Min{l, -^^^r-f—, o'^'fu I / (4-33) 
Suph Sup{b+ |r(7V 

V M-V 

The function (/?_ = > is a subsolution of the original Lichnerowicz 
equation. It is always possible to choose the constant supersolution m large 
enough such that < m on M. 
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2. If we know only that Infb > 0, C/i a subset of M, the problem is with 

Ui 

the supersolutions. We try to conformally transform the metric 7 to a metric 
with positive scalar curvature in Vi, with M — Ui <Z Vi, by considering now 
the equation in Vi and boundary condition, with k some positive number, 

A^e + k9^0, e\9v,^l. (4.34) 

This equation implies that: 

r(7i) = r^(r(7) + /c), if 71 = ^4/(^-2)7. 

with r(7;^) > if 

fc>|r(7)| (4.35) 

The problem 4.34 is of Fredholm type: it has a solution if the homogeneous 
problem, i.e. the equation 4.34 but with boundary condition 9gvi — 0, has 
for unique solution 0, which is the case for A; > 0, if A; is smaller that the first 
eigenvalue. In this case the solution 9 is positive. The conformally trans- 
formed equation (we extend 9 to M, as above) has constant supersolutions, 
all numbers nii such that 

Sup{\r{ji)\ + ai + qi) Sup{ai + qi) 

mi > Max{l, — — , ^] }, (4.36) 

^ Infb /n/^-("+')/("-'V(7i)^ ^ ^ 

A supersohition Lp^ of the original equation, arbitrarily large (but bounded), 
such that Lp_ < is deduced from the inverse conformal transformation. 

a. If r(7) = the condition 4.35 is satisfied by any positive k. 

b. If r(7) = — 1, the condition 4.35 is verified with a number k smaller 
than the first eigenvalue A of the homogeneous problem . It holds that 
A = Cp^ where Cp is the Priedrichs constant of the domain (T^, 7), smallest 
number such that: 

/ \u\^IJ^^<Cf I \Du\^ix^, for all ueViVi). (4.37) 
JVi Jvi 

■ 

Conclusion. Existence and non existence results cover all cases when b, 
that is T, is a constant, since then either Infb > or 6 = 0. A gap remains 

M 

when r is not constant and (M, 7) is in the negative Yamabe class. 
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Remcirk. The method can be apphed in the presence of unsealed 
sources, with practically no change if 6 — ^2 ^ 0. If & — ^2 < there can be 
solutions only if 7 is in the positive Yamabe class (see a construction of sub 
and supersolutions in the case n = 3 in CB 1972). In the general case the 
discussion is more involved, but seems to offer no conceptual difficulty. 

4.4 Scaled sources in W. 

Though discontinuous energy sources are allowed by our hypothesis qi G 
the following theorem, deduced from theorem 8.2, is useful in the presence 
of less regular sources and for the coupling with the momentum constraint 
with low regularity sources. 

Theorem 4.12 // (M, 7) is in the positive Yamabe class the theorem 4-10 
and corollary 4-11 hold with hypothesis on a and qi weakened to a, q\ e W. 

The same result holds for the zero Yamabe class under the additional 
assumption b ^ 0. 

Proof. The problem in the case where a + qi is not bounded above is in 
the construction of a constant supcrsolution m. To construct a non constant 
supersolution, in the positive Yamabe case, we proceed as in Moncricf 1986 
and CB-Moncrief 1994 (case n = 2). For simplicity we write up the physical 
case n = 3. We denote by y?Q = |/o the positive number solution of the 



Such a number exists if r > or r = and h> 0. 

We define one function v € with mean value zero on M, by solving 
the linear equation 



equation 



kvl + Lvl - q^Vo - a = 0. 



(4.38) 



where / denotes the mean value of / on (M, 7) : 




(4.39) 



(4.40) 



The function 



(4.41) 
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is a supersolution because it holds that: 

A^V'^ - /(., = r{ipo - (^+) - a(</?o ^ - (^;^) - gi((^o ^ - ^+^) + Kv^o - '^D 

(4.42) 

hence if r > 

A^^+-/(.,^+)<0. (4.43) 

Since <^_,. > (^q > we can choose a subsolution £ > such that < i < (p^ 
and apply the theorem 8.2. ■ 

When n = 3 our result includes the case p — 2, i.e. 7, G 

5 Coupled system, r =constant, scaled sources. 

When T is a constant the momentum constraint (scaled sources) does not 
depend on (p. We can therefore solve this linear system for X and insert af- 
terwards the result in the semilinear equation for ip. To solve this Lichnerow- 
icz equation we have supposed that 7 G Mf, p > |, and the coefficients 
a, b, qi G except a and gi in the positive Yamabe case. While b and qi 
are given quantities, the function a depends on the given 7 and B, and also 
on the unknown X, solution of the momentum constraint. We leave to the 
reader formulations of variants of the following existence theorem. 

Theorem 5.1 Let 7 G Mf. Let qi and Ji, orthogonal to the conformal 
Kiling fields of {M, 7), be given as well as B & Wf, and r —constant. Then 
the LCBY system with scaled sources admits a solution X,(fi & , (p > 0, 
in the following cases 

1. (M,j) is in the positive Yamabe class, p > |, qi, Ji G If, B G Wf 
and \J\i + qi ^0 on M. 

2. (M,j) is in the zero Yamabe class, p > f , Qi, Ji G L^, r 7^ and 
\J\i + qi ^0 on M. 

3. (M,j) is in the negative Yamabe class p > n, qi & L°° , Ji G LP, 
B G Wf, r ^ 0. 

The solution is unique, except if (M, 7) is in the zero Yamabe class and 
|J|i + |5|+gi + r2 = 0. 

Proof. Under the hypothesis the momentum constraint 3.1 has a solu- 
tion X G VFf , since the vector F given by 3.2, here independent of p and J2, 
is in L^ and is L^ orthogonal to the space of conformal Killing (CK) vector 
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fields of (M, 7). The tensor A = C^^conf^ + -B is then in Wf. The function 
a = A. A is in as soon as p > |. 

1. and 2 The existence theorem 4.12, positive or zero Yamabe class, ap- 
plies if a + (^1 ^ 0. Since 

DiA'^ = 4 (5.1) 
we have A ^ 0, hence a ^ 0, if Ji ^ 0. Therefore \J\i + Qi ^ implies 

3. If p > n the tensor A is in C°, hence the function a is also in C° C L°°. 
We apply the theorem 4.10 and its corollary 4.11. ■ 

6 Solutions with non constant TrK, or (and) 
unsealed sources. 

We will show that there exists a whole neighbourhood of low regularity so- 
lutions with non constant^ r and (or) non zero unsealed sources near such a 
solution with constant r and zero unsealed sources. 

Lemma 6.1 The mapping defined by: 

$ : {x, y) ^ $(x, y) = {n{x, y),M{x, y)), (6.1) 

x={j,T,B,qi,q2,Ji,J2), y={(p,X) (6.2) 

3n — 2 n . ft — 2 9, n+2 

l-i{x,y) = I\^ip-r{^)>~p^a{B,X)ip "-2 -Fgiiy? ^-■^-{02 — r^)(/?»-2 (6.3) 

4n 

(6.4) 

is, ifp > f , a dijferentiable map into the Banach space Bq = x ^ ®U' from 
the open set fl = fli x Q2 of the Banach space B = Bi x B2, defined hy^^ 

fii = Bi n {7 e M|}, = B2 n > 0} (6.5) 

^There are cosmological spacetimes which do not admit surfaces of constant r (Bartnik 
1988) 

^•^To help the reader follow our reasoning we denote by (8) Wf the space of Wf m- 
tensor fields. 
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Bi = {W^ X Wl X (2(g) V^f) X L^' X x^ (gjL^' x^ ^V), (6.6) 

B2 = X C^W^). (6.7) 

Proof. The mapping $ is continuous from Q into LP ^LP because 
</? e W2 C C°, > 0, and DB,dT,a,qi,q2, Ji, J2 E LP as well as and 
r = r(7) when p > |, . It is differentiable on Q because the linear mapping 

S^x,y ■ {Sx, 5y) ^ d<^^^y{6x, dy) = {5n^^y{5x, 6y), 6Mx,y{i!X, dy)) (6.8) 

is a mapping from the Banach space B into the Banach space LP (S>LP for 
each {x,y) e fl. Indeed this linear mapping is given by: 

(37ia;,3/(da;,dy) = A^^(/j-{rH ""'H — ^iT^iV^ — 3i^(o-92)¥'"-4t"^ 



+S'y.D'^ip + S{T{'y)).D(f-6rLp + 6a(p '"-■^ +6qi(p "-^ + {Sq2- Sbjcp^-^ (6.9) 
and 

r«v /r r X /a ri^xi / 2(^—1) 2i+2 2(^+2) n+6 • 

5M^,j,((ya;, 52/) = (A^,eon/5^)* - ( ^ _ V "-^^V + ^ _ ^ V "-^ ■/2)^V^ 



+DjSB'^ + (5(r(7)).S)' - {^^^ — -ip^d'Sr + </^^5J* + (^Ji). (6.10) 

Til 

Embeddings and multiplication properties give as in previous sections the 
announced result. ■ 

Theorem 6.2 Let (g,K), K = B + Cg^conf^ + be a solution of the Ein- 
stein constraints on the compact n-manifold M, with TrK= f a constant, 
g e Mf , X e Wf , B e Wf, P > ^ and scaled sources qi, Ji e LP and no 
unsealed sources. Suppose (M, g) admits no conformal Killing field. There is 
a neighbourhood off,B, qi, q2 = 0, Ji, J2 = in Bi such that the Einstein 
constraints with non constant r and (or) unsealed sources admit a solution 
(g, K) in a x Wf neighbourhood of (g, K). 
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Proof. We have just proved that the mapping $ defined by the left 
hand side of the LCBY equations is a differentiable map into a Banach space 
B2. Saying that the function (p = 1 together with X G W2 is a solution of 
the LCBY system with data 7 = ^ e M|, r a constant, ^1 G L^, , §2 = 0, 
B G Wl, Ji E If, J2 — is to say that, with the notations of the previous 
lemma: 

$(x,y) = (6.11) 

.By the implicit function theorem^^ there exists of a solution y ={if, X) G ^^2, 
i.e. if eW^, ^ > 0, X e Wi, when x = e Mf , r G Wf, B G V^f, gi, gs, 
Ji, J2 G L^} is in a sufficiently small neighbourhood in the Banach space B 
of X, under the condition that the partial derivative of $ with respect to y 
at (x, y) is an isomorphism of Banach spaces, that is if the linear system 

{^-g^conf^Xy (6.12) 

and (recall that = 1) 

AC- r- 3?7.-2 n n + 2-, . , 

A,5^ - + ^« + + = ^ (6.13) 

has one and only one solution 6X G^ ®W^2 ' ^ ^2 ^ 
ijj E LP. This result holds for 6.11 if g admits no conformal Killing field. 
It holds for 6.12 if the coefficient of Sep, which belongs to If, p > | by 
hypothesis, is positive or zero, and not identically zero. It is obviously so 
when f > 0. To treat general values of f we remark, as O'Murchada and 
York 1974 in their study of linearisation stability, that since {g, K) satisfies 
the hamiltonian constraint, hence the Lichnerowicz equation with ^ = 1, it 
holds that: 

-f + a + gi-6 = 0. (6.14) 

The linear equation 6.12 reads therefore 

2(n — 1) 4 - 
A-g5^ - {2a + ^ -^qi + -h}5p^ = ^ (6.15) 

The coefficient of the lower order term 5ip is positive or zero, and not iden- 
tically zero, if a, gi, and h are not all identically zero. ■ 

The result does not extend in a straightforward way^^ to the case where 
g admits conformal Killing fields because M. does not takes its values in the 
space of vectorfields orthogonal to such fields. 

^^See for instance CB and DeWitt I p. 91. 

^^An analogous problem occurs in linearization stability (see Moncrief 1975) 
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7 Appendix A. Second order linear elliptic 
systems. 

7.1 Linear systems 

Let (M, e) be a smooth compact orientable riemannian manifold. Denote by 
d the covariant derivative in the metric e. A second order Uncar differential 
operator from sections w of a tensor bundle E over (M, e) into sections of 
another such bundle F reads 

2 

Lu = ^akd^u (7.1) 

fc=0 

with ttfc a linear map from tensor fields to tensor fields given also by tensor 
fields over M. In local coordinates it reads: 

{Lu)^ = 4;'Jd^y + 4%diu'' + a^gu^ (7.2) 

The principal symbol of the operator L at a point a; e M, for a covector 
^ at X, is the linear map from to determined by the contraction of 02 
with represented in coordinates by the matrix 

= <4;iiu, (7.3) 

The operator is said to be elliptic if for each x E M and ^ G T*A4 its 
principal symbol is an isomorphism from E^ onto F^ for all ^ 7^ 0, i.e. the 
determinant of the matrix does not vanish. If there is a covering of M 
such that this determinant is uniformly bounded away from zero the system 
is said to be uniformly elliptic and the bound is called the ellipticity constant. 
The Laplace operator in a riemannian metric 7 on M acting on scalar 
functions u has principal symbol Y^CiCj: it is elliptic. 

The Sobolev spaces W'^ on (M, e) are defined as closures of spaces of 
smooth tensor fields in the norm 

\\f\\wf^{[ E l^'/IVe}^ (7.4) 

where d, \.\ and fi^ denote the covariant derivative and the pointwise norm 
in the metric e. 
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We denote by C any positive number depending only on (M, e) : volume, 
Sobolev constant... 

A metric 7 is said to belong to if it is properly riemannian and 7 G W^. 
We will always suppose that a > -, then 7 e is continuous on M and 
is an open subset of W^. The contrevariant associate 7* is also continuous 
and the volume elements /i^ and //g are equivalent. 

We recall the following theorem. 

Theorem 7.1 (Doughs and Nirenberg) Let fl be a bounded open set of K", 
let L be an homogeneous elliptic operator: 

Lu = amd'^u (7.5) 

where d is the usual partial derivative and a^n is continuous and bounded on 
then the following estimate holds, for any p > 1: 

\\u\\w^ < Cam{\\Lu\\ip + \\u\\ip} (7.6) 

where Wf denotes Sobolev spaces defined by the euclidean metric on fi. The 
number Ca^ depends only on Q the ellipticity constant, the norm of 
and its modulus of continuity. 

We will deduce from this theorem the following one^^. 

Theorem 7.2 Let (M, e) be a smooth compact orientable Riemannian man- 
ifold. Let 

2 

Lu=J2 (^kd'^u (7.7) 

fe=0 

be a second order elliptic operator on (M, e). Suppose the coefficients of L 
are such that: 

a2 e Wi, ai e Wf, ao e L^, p> -. 

Then 

1. The operator L is a continuous mapping Wf ~^ ■ 

^■^Lowering the regularity of the coefficients of the theorem proved in CB-Ch, it can be 
extended to higher order elUptic operators and possibly lower values of p for the solution, 
but wc will need the restriction p>n/2 when working with the non linear Lichnerowicz 
equation and treat here only this case for simplicity. 
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2. The following estimate holds: 

IIwIIh/p < Cl{\\Lu\\lp + II^^IIli}. (7.8) 

The number Cl depends only on the norms of the o^s in their respective 
spaces and the ellipticity constant 0/02- 

Corollciry 7.3 // morevover 

t/ien L is a continuous mapping VFf_,_2 — >^ and for all u e W^f+2 /loMs 

Mwl^^ < Cl{\\Lu\\w^ + \\u\\li}. (7.9) 

with Cl a number depending only on the norms of the a^s and the ellipticity 
constant of 02. 

Proof. By the Sobolev embedding theorem we know that 02 belongs to 
the H51der space C°'°^ since C C^'" if p > f . Therefore aa^^u e if 

If ai G VFf then aidu G due to the continuous multiphcation Wf x 
Wf LP when p > | and it holds that 

\\aidu\\LP < C\\ai\\wp\\du\\wp (7-10) 

If ao G I/^ then cqM G whatever be m since u G C°'". 

To prove the inequality 7.8 using the Doughs - Nirenberg theorem we 
proceed in two steps. 

a. We treat the principal part by considering a covering of M by a 
finite number of charts with bounded domains Qj, and a partition of unity, 
functions < 0^ < 1 with compact support in fit such that ^^^j = 1. We 
have u = "^Ui, Ui = u(f>^. In any of the charts e is uniformly equivalent to 
the euclidean metric; we denote by d the derivative in the euclidean metric 
(i.e. the usual partial derivative), by the associated Sobolev norm. For 
fields with support in a chart the norms Wf and Wf are equivalent, they are 
uniformy equivalent for the various charts, since there is a finite number of 
them. The Doughs - Nirenberg theorem gives then the existence of a constant 
Cj, depending only on fij and on a2 through its C°'" norm and the ellipticity 
constant of its representatives, such that 

\\Ui\\w^ < Ci^a2{\\LUi\\LP + \ \u\\lp}, LUi = a25^Wj. (7.11) 
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The subadditivity of norms together with this inequahty imply that, with 
the maximum of the C'^s, 

\\u\\wi < X] < ^{||Lttj||LP + (7.12) 

i i 

On the other hand, we have an identity of the form, with S{i) smooth coef- 
ficients hnear in the connection of e : 

Lui = a2d^Ui = a2{d^Ui + S{i)dui} (7-13) 

Using the Leibniz rule we obtain that 

dui — udcf)^ + du(l)i, d'^Ui — d'^u(f)j^ + 2dud(f)j^ + ud'^cf)^. (7-14) 

Since 0< 0j < 1 we deduce from 7.14 that 

\a2d'^Ui\ < \a2d^u\ + \a2{2dud(t)i + m^VJI- (7.15) 

Combining previous inequalities we find that there exists a number C de- 
pending only on (M, e), such that 

||Ltii||LP < C{\\a2d'^u\\LP + ||«2||co||«||wf}- (7-16) 
The inequality 7.12 implies then that: 

\\u\\w^ < Ca2{\\a2d'^u\\LP + ||m||lp}. (7.17) 

b. We now consider the general linear operator L, with coefficients satis- 
fying the given hypotheses: 

Lu = a2d'^u + aidu + QqU. (7.18) 

Since 02 G ' with p > |, it belongs also to a Holder space C°'", < a < 
2 — ^, and it holds that, with C a Sobolev constant of {M, e), 

||a2||co,c < C||a2||uA|'. (7.19) 
The inequality 7.17 implies that 

||it||w|' < C'o2{||-^'*^||lp + \ \aidu + aou\\LP + \ \u\\wp}- (7.20) 
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To estimate ||ai9w||iP we use the Holder inequality 

\\aidu\\LP <\\ai\\LPi\\du\\L<,, — + - = -. (7.21) 

Pi Q P 

The Sobolev embedding theorem gives that 

Tip 

\\ai\\LPi < C\\ai\\wr, Pi = • (7.22) 

^ n — p 

Taking this value of pi gives that 

11 11 

q P Pi n' 

and 

\\aidu\\LP < C\\ai\\wf\\du\\Ln. (7.24) 

On the other hand: 

IIoo^IIlp < ||oo||lp||w||l°° ('''■^s) 

By the Gagliardo-Nirenberg interpolation inequality (see for instance 
Aubin p. 94, CB-DM II p.385) it holds that 



(7.23) 



\\&>u\\l.<C{\\u\\'^,\\u\\Y,^ (7.26) 

with: 



A=^{^ + ---} (7.27) 
2 n p q 

if 

< A < 1. (7.28) 
Applying the equality 7.27 to the case j — l,q — n gives 

< A = — < 1 since p> - (7.29) 

2p 2 ^ ' 

Therefore, with this value of A, it holds that: 

\\du\\Ln<C\\d\\\l,\\u\\^^,^. (7.30) 
The same procedure applied to the case j — l,p — q^\— \ gives that 

\\du\\Lv < C\\d'^u\\l^\\u\\l^. (7.31) 
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The interpolation inequality 7.26 with q — oo and j — gives, with again 

2p 

||ti||Loc < C\\d'^u\\^Lp\\u\\l-^. (7.32) 
For any £ > 0, by elementary calculus, it holds that 

WdMliAM'LP^ < ^£\\dM\Lp + ^^_^j^;,/,_J |^IUp (7.33) 

hence there exists a number C such that 

\\aidu + aou\\LP < C{\\ai\\wP + \ \ao\\Lp}{s\\u\\w^ + C^\\u\\lp} (7.34) 

Using this inequality together with 7.20 we see that we can choose £ > 
small enough, depending on the bound of ||ai|||^p + HooUlp and Ca^ such 
that the inequality 

\\u\\w^<Cl{\\Lu\\lp + \\u\\lp} (7.35) 
is satisfied. But by the definition of integrals and norms it holds that 

ii^iIlp < Wuwt^'^^^wuw'jir < (7.36) 

Therefore if > |, for any e > there exists e and such that: 

\\u\\lp < I l^illj^p^^^^l 1^*11^1^ < ^II^IIm/I + QIImIIli, (7.37) 

we can again choose e so that the announced inequality holds. ■ 

Proof, of corollary. This regularity statement can be proved by stan- 
dard recursive arguments. ■ 

Theorem 7.4 Under the hypothesis of the theorem 7.2 for s = 0, and its 
corollary 7.3 in the case s > 0, it holds that: 

1. The operator L maps W2_^_g into with finite dimensional kernel and 
closed range. 

2. If L is injective on then there is a number Cl such that for each 
u in Wgj^2 the following inequality holds: 

\Mwi^^ < Cl\\Lu\\wp . (7.38) 

3. If the formal adjoint *L of L satisfies the same hypothesis as L and is 
injective, then L is surjective from Wf^^ onto W^, hence an isomorphism if 
also injective. 

Proof. It follows usual hues. ■ 
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7.2 Isomorphism theorem for the Poisson operator. 

The Poisson operator in a riemannian metric 7 on {M, e) acts on scalar 
functions u and reads, with S^j the difference of the connections of e and 7 : 

A^u -au = f^id'^jU + S^jdhu) - au. (7.39) 

We suppose that 7 e M|, p > f , then 7 e C°'". We have denoted by Mf 

the open subset of symmetric 2-tensors in VFf which arc properly rieman- 
nian metrics, namely such that they admit a positive ellipticity constant, 
/n/det(7jj) > 0, where 7^^- are the components of 7 in a finite number of 

charts M/ covering M. The coefficient 0,2 = 7*, the contrevariant tensor asso- 
ciated to 7, is then in C°'" with a positive ellipticity constant. The coefficient 
ai = Y^Sij is in Wf because d^d^ e i/ if 7 G P > f • The hypothesis of 
the theorem 7.2 are satisfied if a e i/. Hence: 

Lemma 7.5 If j E Mf, p > f , i/ie Poisson operator is a continuous map- 
ping H/f ^ LP if a e LP. 

We now prove a uniqueness lemma^^. 

Lemma 7.6 Let 7 e Mf , p > | , a G L^, ^/ie Poisson operator — a is 
injective on if a > 0, a ^ (i.e. a > on a subset of M of positive 
measure). 

Proof. The proof for functions results from the maximum principle 
if the coefficients are bounded. In the more general case one proceeds as 
follows: on a compact manifold, if m„ is a function and 7^ G C^, the 
following identity is obtained by a straightforward integration by parts 

/ Un{A^^Un)n^^ = - {j'^diUndjUn)fly^ (7.40) 

Jm J m 

This identity applied to a solution of tt„ — aw„ = 0, implies that 
diUn = on M, and = on the subset a > 0. Hence rt„ ^constant, and 
= . 

Suppose now that 7 G Mf, p > f , and u G VKf . Since is dense in W2 
we can approximate 7 and m, and a by smooth sequences 7„, m„ converging 

^"'The conditions p>n/2 for 7 and for a are not necessary, but they are sufficient. 
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to 7 and respectively in Mf , VFl. We already know by the lemma 7.5 that 
A^^M„ converges in L'^ to AyU. Since the measures and are equivalent 
to /Xg and the continuous embedding x always holds as well 

as Wf X Wf C -L^ if p > |, the integrals on both sides of 7.30 converge to 
limits, and it holds that: 

/ u{A^u — au)fj,^ = — / {Y-'diudjU + au^)ijL^ (7-41) 

from which follows it = if A^u — au — and a > 0, a ^ 0. ■ 

Theorem 7.7 The Poisson operator — a on scalar functions in a metric 
7 on a smooth compact riemannian manifold {M,e), with 7 e Mf, p > f , 

a G L^, is an isomorphism from W2 onto U' if a. > 0, a ^ 0. 
It is an isomorphism Wg_^2 ~^ > ^/ addition a e VFf . 

Proof. Under the given hypothesis the Poisson operator is selfadjoint 
relatively to the metric 7 : two integrations by part for smooth enough 7, 
u and V, and taking limits, show that under the given hypothesis, as in the 
proof of the lemma 7.6, the following identity holds for n, f e Wf : 

/ v{A^u — au)iJ,^ = / u{AjV — av)iJ,^. (7.42) 

The isomorphism theorem 7.4 applies therefore to our Poisson operator. ■ 



7.3 Generalized maximum principle. 

One knows by the classical maximum principle^^ that if a solution u & oi 
the inequality with bounded coefficients and a > 

AjU — au < 0, [respectively A^u — au >0] 

attains a minimum A < [respectively a maximum A > 0] at a point of M 
then u = X on^^ M. Therefore A^n — an < on M implies n > on M if 
a > on M. We will need to obtain low regularity solutions the following 
generalization of the maximum principle. 

^'"'Protter and Weinberger 1967, p. 64. 
Also if S is a bounded domain of M with smooth boundary a minimum < [respectively 
a maximum > 0]must be attained on the boundary. 
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Lemma 7.8 IfuE ; P> 2 ^(^tisfies the equation 

A^u -au = -f (7.43) 

with 7 e M|, a e L^, a > k, with k > a number , and f > , then u >0 
on M. 

Proof. The lemma holds by the classical maximum principle if u G C^, 
7 e M| a e C°. We approximate 7 G M| by 7^ G Mf , a by a„ G C°, a„ > A; 
and / G i7 by /„ G C*^, /n > 0. We know by the isomorphism theorem that 
each equation 

Ln = Ay^Uji Cbn^n — fn 

has one solution w„ G Wl C C^, and -u^ > 0. We know that there exists a 
number Cl depending only on the W2 norm of 7„, hence of 7, and the 
norm of a„, hence of a, such that: 

llMnllw^l' <CL{\\fn\\LP+\\Un\W}- (7.44) 

We also know that, since L„ = A^^ — a„ is injective on W^, there exists a 
number Cl„ such that: 

<ClJ|/„||l., (7.45) 

however the number Cl„ does not depend only on the ^ norms of 7„ 

and a„, hence is not a priori independent of n. In order to obtain our positivity 
result by a continuity argument we must estimate C^^ independently of n. 
We proceed as follows. 

The equality 7.41 shows that a solution u G p > |, of 7.43 satisfies 
the integral equality 

- / u{A^u - au)fj,^ = / Y^diudjU + au^)/!^ = / uf/i^. 
J Jm J m 

which implies if a > A; > 

/ u'^ii^ < k-' [ ufii^ (7.46) 
Jm Jm 

This inequality and the uniform equivalence through a constant C-y (depend- 
ing only on the norm of 7 and its ellipticity constant) of the quadratic 
forms 7 and e, together with the Holder inequality give: 

M\l2<C^,kM\Li\\f\\L,', ^ + ^ = 1' C^,k = C^k-' (7.47) 
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We have supposed f & LP. 

Suppose first that p > 2. We have then, since M is compact, f & L'^ with 

II/IIl^<C^II/IIl. (7.48) 
We choose then q' = 2, q = 2 and we deduce from 7.46 

\\u\\l^<C\\u\\l^<C^Mlp. (7.49) 

Hence for a solution it of a Poisson equation with coefficients 7 e Mf , a E LP., 
p>|,a>A;>Oan inequahty 

\Mwi'<CL,k\\f\\Lr', (7.50) 

where CL,k depends only on the Wf norms of 7 and a, and the value 

oik. 

Suppose that p <2 (inequahty compatible with p > f only if n < 4) . We 
then take q' — p, but we have q > 2, since q"^ — 1 —p~^. We use the equahty 

\\u\\L.<C\\u\\%f/''\\u\\'^,' (7.51) 

to obtain that 

IHt'<C||«lli^|'^^''ll/lk^ (7.52) 

hence 

II^IU. < < Ce\\u\\^. + CM\l^ (7.53) 

from which follows again an inequality of the type 7.50. 

We apply a similar inequality to the difi'crcnce Un — Un-i which satisfies 
the equation: 

^■yJUn - Un-l) - a^iu^ - Un-l) = (7.54) 

with 

Fn = (A^„ - ^7„_i)Mn-l " (/n " fn-l) + (^n " a^-ljUn-l. 

The expression of as linear in the second and first derivatives of u, the 
usual embedding and multiplication properties, and an inequality of the type 
7.50 satisfied by w„_i, show that F^ G with an L^ bound of the form 

I WuWlP < CjJ,k{\\ln - In-lWwS + Wfn - Iu-iUlp + IK - a„_i| \lp}- (7.55) 
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where C^,/,fc denotes any positive number depending only on k, the W2 norm 
of 7. and the norms of a and /. The inequahty 7.50 apphed to 7.54 gives 
therefore an inequahty of the form 

\\Un - Un-l\\w^ < Cj,f,k{\\lfn " 7n-lllw| + ll/n " /n-lHif + " On-lHLf}- 

We deduce from this inequahty that the sequence it„ converges in W2, hence 
in C°, to a function u G W2 when converges to / and 7„ to 7. This 
function u is the unique solution of the Poisson equation, and since Un > 0, 
the same is true of -u^''. ■ 

8 Appendix B. Equation A^cp = /(x, (p). 

We consider equations of the form, with (p a scalar function 

A,(^ = /(.,^). (8.1) 

where / : {x,y) 1— > f{x,y) is a function on M x /, with / an interval of R, 
smooth in y. To be specific, and in view of application to the Lichnerowicz 
equation, we suppose that / is a finite sum: 

f{x,y)^ J2 ai{x)y''', (8.2) 

1=1,.. .N 

where the exponents Pj are given real numbers. Hence / is a smooth function 
of ?/ if 7/ > £ > 0. 

It is to solve this semi- linear equation that we have supposed p > The 
space W2 is then an algebra, is in W2 C C° for any P if it is so of ip and 
if ip is positive. 

To solve 8.1 we apply a method^^ of successive iterations obtained by 
resolution of linear problems, whose bounds and convergence are obtained 
through the use of sub and super solutions. 

^^Recall that it is to obtain a family of smooth functions u„ approaching u that we had 
to introduce the positive k. Generahzed maximum principles can probably be obtained 
by other methods. 

^^The solution of equations of the indicated type on compact manifolds was first ob- 
tained, in Holder spaces, by using the Leray-Shauder degree theory, together with bounds 
by sub and super solutions in CB and Leray 1972. It was applied to the constraints by 
CB 1972 (imscaled sources), O'Murchada and York 1974 (scaled sources), extended to 
Hg spaces in CB 1975. The iteration method (Kazdan and Warner 1985) was applied to 
the constraints by Isenberg 1987. In the following we weaken the regularity hypothesis 
previously made. 
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definition 8.1 1V| function ip_ is called a subsolution of Aj(p — /(., cp) 
if it is such that on}^ M 

A^ip_>f{.,^_). (8.3) 
An W2 function t/?^ is called a supersolution if 

A^V9+</(x,(^+). (8.4) 

Theorem 8.2 (existence) The equation A^Lp = /(., (p) on the compact man- 
ifold (M,'~f) admits a solution ip G W2, p > f , if-' 

a. 7 e M|, p > f , ai e LP, 1=1,. ..,N. 

b. The equation admits a subsolution ip_ and a supersolution (p_^, both in 
W2, < £ < (p_ < (p_^_ < m. The solution is such that (p_ < (p < (p^- 

Proof. The proof follows the same lines as the standard proof, but 
uses the generalized maximum principle (lemma 7.8). We define successive 
iterates by solution of the elliptic equation, linear in (p^ when fn-i is known, 

^I'Pn - a'Pn = f{; ^n-l) " a'Pn-l, (8-5) 

with a e LP, a > k, k a. strictly positive number to be determined later, 
introduced now to allow the application of the theorem 7.7. 

We take (p^ = (fi_ to start the iteration, i.e. we define (pi. by solving the 
equation 

A^(/Ji - a(/?i = /(.,V?_) - av?_, (8.6) 

due to the hypothesis made on j,a,f and (p_, this solution 99^ is in lyf- 
Moreover we have by the definitions of (p^ and (p_ 

A^iip, - ip_) - a{^^ -ip_)<0 (8.7) 

The conclusion (p^ — (p_ > on M follows from the generalized maximum 
principle. 

We now use the definition of (p_^. to obtain: 

A^(V7+ - ip^) - a((/7+ - v?i) < f{x, v?+) - f{x, ip_) - a(v?+ - ip_) (8.8) 

which gives using the mean value theorem (CB-DM I p. 78) 
^^Inequality satisfied for functions in L^, i.e. almost everywhere 
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A^(^+-^i)-a((/?+-¥?i) < [ {f^{x,ip_+t{^+-^_))-a}dt (8.9) 

It is sufficient to choose the function a such that 

a = Supi<y<^\f;^{.,y)\ + k, k>0 (8.10) 

to have on M that a > A; > and 

A^((^+ - ip,) - a(^+ - (^i) = F < 0, (8.11) 

since on M it holds that <£<</?_ + 1(</7^ — < m. For example if the 
function / is given by 8.2, then: 

i/;(.,2/)i < J2 \^iPi\y'''~' < E WuPiM'^'^~'+J2 Wi-Pi-\^'''-~'+\aio 

1=1,.. .N 1+ 7_ 

(8.12) 

with P/_|_ > 1, Pj_ < 1, = 1. Under the hypothesis made on a/ the 
maximum principle holds, hence Lpi<(p_^_on M. 

Suppose that there exists (/j^ G VFf, m — 1, ...,n — 1, such that on M 

^-<^m< ^+ (8.13) 

then (p^ exists on M, solution of 8.5, also in . Suppose moreover that 

V^n-l - ¥>n-2 > (8.14) 

Then the equations satisfied by y?^ and (p^-i imply 

^liVn - Vn-l) - (^iVn " "Pn-l) = i^n - a)(^n-l " fn-2) (8-15) 

where 6„ is the function in 

bn= f!p{-: + t{^n-l " <^n-2)dt. (8.16) 

JO 

It results from the definition of a that: 

hn-a<-k<Q (8.17) 

this inequality together with v^n-i ^ ¥'n-2 insures that the right hand side 
of (8.15) is negative, hence (/?_ < ^fin-i — Vn- To show that the sequence is 
bounded above by (/?^ wc write again an inequality of the form 

A^((^+ - if^) - a{ip^ - ipj < f{ip+) - /(^„_i) - a((^+ - ^„_i) (8.18) 
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from which the conclusion follows. We have proved the existence of the 
sequence G Wf) with <(/?_< < hence Hv^nllip uniformly 
bounded by Due to the defining equation 8.5 llv^nll^yl ^1^° uni- 

formly bounded. Since functions in W2 are equicontinuous we can extract 
(Ascoli - Arzela theorem) from the sequence a subsequence, ip^, which 
converges in C° norm, to a function (p e On the other hand the sequence 

of continuous and positive functions is increasing and bounded above by 
(p^, hence it is pointwisc convergent to a funtion (/?, which coincides with 
the previously obtained (p. The whole sequence (not only a subsequence) is 
therefore convergent. 

We show that the limit is in fact in by using the elliptic estimate 
applied to the equation 8.15 

II fn - lk|< C^A\{bn - a)(^n-l " \\lv (8.19) 

which implies: 

II Vn - Vn-1 \\w^< C7,all(^n " o)\\Lp\Wn-l " 'Pn-2) lico, (8-20) 

The sequence tp^ converges in W2 (to ^p by uniqueness of limits) since it 
converges in C*^, and the hn — a are uniformly bounded in LP norm. The 
limit satisfies the equation in L'^ sense. 

The solution we have constructed depends on the choice of the initial (p^. 
We have taken y?_ as an initial (p^. We could have taken (pj^_ and construct 
a decreasing and bounded below sequence which converges to a limit ^. 
We cannot in the procedure used start from an arbitrary <^q, with p)_ < 

^ <^+) because we will not know if (p^ satisfies the same inequalities. 
A uniqueness theorem is easy to prove in the case where y 1-^ fi-,y) is 
monotonically increasing. A more general uniqueness theorem (4.4) holds 
for the Lichnerowicz equation. ■ 

Aknowledgement. I thank the referee for useful criticisisms. 
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